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The authors of the preceding Comment [G. M. D'Ariano 
and C, Macchiave llo Phys. Rev. A (preceding comment), 
quant-ph/9701009] tried to reestablish a 0.5 efficiency bound 
for loss compensation in optical homodyne tomography. In 
our reply we demonstrate that neither does such a rigorous 
bound exist nor is the bound required for ruling out the state 
reconstruction of an individual system [G. M. D'Ariano and 
H. P. Yuen, Phys. Rev. Lett. 76, 2832 (1996)]. 
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There is little doubt that the compensation of detec- 
tion losses is a numerically delicate procedure that is ex- 
tremely sensitive to experimental inaccuracies. We have 
shown UH, however, that no clear efficiency bound exists 
beyond that loss compensation is impossible, in contrast 
to a statement in an earlier paper p[. In the preceding 
Comment [Q D'Ariano and Macchiavello tried to reestab- 
lish a 0.5 bound for the overall efficiency r\. 

In this Reply we show that their analysis is incomplete 
and that still no in-principle bound exists. Furthermore 
we point out that the existence of such a compensation 
bound does not follow directly from the arguments given 
in Ref. || where the impossibility of measuring the state 
of an individual quantum system was proven. 

What is the problem? We have shown jl| that a gen- 
eralized Bernoulli transformation describes the influence 
of detection losses on the density matrix {m\p\n) that is 
reconstructed from measured data obtained e.g. from op- 
tical homodyne tomography ^| . We have treated the loss 
process separately from the particular detection scheme. 
The inversion of the Bernoulli transformation [Q pro- 
duces the unperturbed density matrix of the signal 
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We have shown Q that the convergence of the series ([!]) 
is not necessarily restricted to the range 0.5 < r\ < 1, al- 
though the matrix elements Bj^^ 1 ) of the inverse trans- 
formation are divergent if 77 < 0.5 0. Furthermore ||] 



we could employ multiple runs of compensation or the 
analytic continuation of the series ([!]). In this case no 
convergence bound on 77 does exist and hence loss com- 
pensation is possible in principle. 

In practice, of course, experimental imperfections will 
affect the compensation of detection losses || . We have 
never claimed that loss compensation is easy ||, in fact, 
we have stressed Q that for rj < 0.5 other errors are 
amplified, e.g. the effect of any uncertainty in rj itself. 
D'Ariano and Macchiavello Q considered the influence of 
statistical errors for homodyne measurements. They as- 
sumed a finite number N of experiments, i.e. a finite sta- 
tistical ensemble of individual quantum systems. In this 
case the reconstructed density matrix (m\p meas (N)\n) is 
an estimation of the matrix (m\p meas \n) with statistical 
error bars calculated according to Ref. jtO|. However, 
as a fundamental axiom of quantum mechanics, the es- 
timation must tend to the ensemble average when the 
ensemble size approaches infinity. Therefore, if we keep 
the cut-off jm in the series ([[]) at an arbitrarily large, 
but fixed value and increase the number N of experimen- 
tal runs, we must approach the correct result {m\p S i g \n), 
with an arbitrarily small, fixed systematic error that de- 
pends only on the cut-off. D'Ariano and Macchiavello Q] 
did exactly the opposite. They fixed the size A?" of the sta- 
tistical ensemble and increased the cut-off ju , and found 
that the series diverges for rj < 0.5. Our Fig 1. illustrates 
the influence of varying both j'm and N. We see clearly 
that the order of the limits jm 00 and A^ — > 00 is im- 
portant. We also see that quite a large number A^ of runs 
is required to produce faithful data for compensating a 
low efficiency 77. D'Ariano and Macchiavello [Q discussed 
certainly an interesting aspect of loss compensation but, 
as we have seen, their analysis was incomplete. 

Is the compensation bound |^] relevant for some more 
fundamental features of quantum mechanics than merely 
technical points of measurement technology? Does the 
50% bound (D ru le ou t the state measurement of an in- 
dividual system (|]? If the bound existed it clearly would, 
as was pointed out in Ref. Q . From the conclusion, how- 
ever, does not follow the premise. The impossibility of 
measuring the wave function does not imply that the 50% 
compensation bound exists. 

The problem discussed in Ref. (5| is again a matter 
of performing limits in the right order. Suppose one 
taps a series of A^ probe beams from an individual light 
mode and performs a state reconstruction using the A^ 
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probes as a statistical ensemble. The effect of tapping, 
i.e. beam-splitting, is equivalent to detection losses |Ly] 
with an efficiency i] scaling like N^ 1 . Therefore, when 
we attempt to reconstruct the quantum state of the in- 
dividual light mode we should employ infinitely many 
probes, yet with infinitely poor efficiency. Not surpris- 
ingly, we cannot compensate the losses in this situation. 
A general 50% efficiency bound is much too much to be 
required for such a delicate matter. 

In conclusion, compensation for low overall detection 
efficiency is numerically difficult. The value 0.5 for r\ 
plays clearly a crucial role || because at this value the 
matrix elements in the inverse Bernoulli transformation 
become unbounded Our analysis shows, however, 
that 0.5 is neither a rigorous bound for compensating 
losses in optical homodyne tomography, nor is this bound 
required for ruling out the state measurement of an indi- 
vidual system ||. 
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FIGURE CAPTION 

Plot of the loss-compensated density-matrix element 
Poo with varying cut-off Jm and ensemble size N. We 
employed the same thermal state with n = 2 as in the 
preceding Comment j| and used an efficiency rj of 0.48. 
We performed Monte-Carlo simulations to model a re- 
alistic experimental situation. First, we reconstructed 
the density matrix (m\p meas \n) from N runs of the com- 
puter experiment using the method developed in Ref. 
JlOj l . Then we performed the loss compensation (|l|) with 
varying cut-off jm ■ We found that for a given jm the re- 
constructed matrix elements poo do approach the actual 
value of (0|p s j s |0) = 0.33 for increasing numbers N of 
runs, apart from a small systematic error. On the other 
hand, if we keep the number N of runs constant and in- 
crease the cut-off 2m the matrix element diverges JiJ , as 
can be seen from the behavior of poo for N = 10 3 . Thus, 
the order of the limits jm — * oo and N — > oo is vital to 
the loss-compensation procedures [QJ|] . 
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